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Stacking interactions in single stranded nucleic acids give rise to configurations of an annealed 
rod-coil multiblock copolymer. Theoretical analysis identifies the resulting signatures for long ho- 
mopolynucleotides: A non monotonous dependence of size on temperature, corresponding effects on 
cyclization and a plateau in the extension force law. Explicit numerical results for poly(dA) and 
poly(rU) are presented. 

PACS numbers: 87.15.-v, 61.25.Hq, 87.14.Gg 
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Single stranded nucleic acids (ssNA) experience stack- 
ing interactions Q]. These favor parallel orientation of 
adjacent aromatic rings of the bases giving rise to rigid 
helical domains. Thus far, the possible coupling of stack- 
ing and the elasticity of ssNA received little attention and 
its existence became recently a subject of debate 0,130 
The issue is further complicated because the relevant 
thermodynamic and structural parameters reported vary 
widely. In this letter we present a theoretical analysis of 
the configurations and elasticity of ssNA subject to stack- 
ing and identify qualitative effects signalling the coupling 
of stacking with the chains' elasticity. Clear signatures 
of stacking are discernable in long homopolynucleotides, 
under high salt conditions when loops do not form and 
electrostatic interactions are negligible. There are two 
primary effects. One is the occurrence of a minimum 
in the radius of the chain, R, as the temperature, T, is 
varied. This leads to corresponding effects on the cy- 
clyzation of the chains. The second effect is a plateau in 
the extension force law of the ssNA subject to tension /. 
The analysis utilizes a model for the helix-coil transition 
in helicogenic polypeptides [j| modified to allow for the 
weak cooperativity of stacking. It focuses on the differ- 
ences between ssNA that stack strongly, polydeoxyadeny- 
late (poly(dA)) and polyriboadenylate (poly(rA)), and 
those that exhibit weak or no stacking, polyribouridy- 
late (poly(rll)) and polydeoxythymidylate (poly(dT)). 
Such long homonucleotides, comprising thousands of nu- 
cleotides, can be readily synthesized enzymatically 0. 
The results suggest that under physiological conditions, 
these effects are important only for poly(dA), poly(rA) 
and ssNA containing large A domains. In the absence 
of such domains, stacking effects become noticeable at T 
lower than 10°C. 

The extension force laws of A ssDNA, as measured in 
optical tweezers or AFM experiments, do not reveal sig- 
natures of stacking Q • The results can be rationalized by 
considering ssDNA as freely-jointed chain characterized 
by a single Kuhn length. This basic picture is augmented 
to allow for loop formation |8| and for electrostatic inter- 
actions 0- However there is no evidence for large A 
domains in A ssDNA and these measurements were car- 
ried out in ambient T thus precluding, as we shall dis- 



cuss, significant stacking effects. Stacking effects were 
reported @,Si 

in "molecular beacons" 101. These are 
short ssDNA chains capable of forming stem-loop struc- 
tures. One end carries a fluorophore and the other a 
quencher. Accordingly an open hairpin fluoresces and 
a closed hairpin is quenched. The fluorescence inten- 
sity and its fluctuations allow to extract the fraction of 
hairpins and the opening and closing rate constants. Ex- 
periments by the Libchaber group revealed differences in 
cyclization behavior of poly(dT) and poly(dA) loops that 
were attributed to stacking and its effects on the rigidity 
of the chains 0, El] ■ This interpretation was disputed by 
Ansari et al. who ascribed the effects to transient trap- 
ping of misfolded loops while arguing that both poly(dT) 
and poly(dA) behave as flexible polymers Our anal- 
ysis does not pertain directly to this system since it con- 
cerns long chains. However, the resulting predictions 
identify clear signatures of stacking when misfolding is 
not an option. 

Stacking involves interactions between nearest neigh- 
bors and is thus non-cooperative or weakly coopera- 
tive pj. It involves a broad transition between the 
stacked, helical state obtained at low T, and random- 
coil configurations at high T. At intermediate T, ssNAs 
comprises of stacked domains interspaced with "melted" , 
random-coil ones. The polydispersed domains undergo 
dynamic equilibrium and the overall behavior is of an- 
nealed rod-coil multiblock copolymer. The strength 
of the stacking interactions vary with the identity of 
the bases. It is strongest between adenosines (A) and 
it is weakest among uridines (U) and thymines (T). 
The interactions between chemically different bases are 
weak. Thus, stacking is most pronounced in poly(dA) 
and poly(rA) while it is weak in poly(rll) and in het- 
cropolynucleotides without extended A domains. There 
is evidence that stacking does not occur in poly(dT). 
In every case, stacking is insensitive to the concentra- 
tion of salt. We consider ssNA homoplymers comprising 
N 3> 1 identical monomers, nucleotides, of which ON are 
stacked. The stacked bases form yN helical domains. In 
comparison to the non-stacked bases, the excess free en- 
ergy of each stacked monomer is A/ = Ah — TAs. Ah 
reflects the enthalpy gain associated with the stacking 
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while As allows for the loss of configurational entropy 
due to the parallel orientation of the stacked bases. The 
reported values of Ah and As vary widely. For poly(rA) 
Ah ranges between —3 to —10 kcal/mole while As val- 
ues span the range —10 to —27 e.u.. In our calculations 
we will use two sets of values: Ah = —13 kcal/mole 
and As = —40 e.u. as reported for poly(dA) as well as 
Ah = —2.7 kcal/mole and As = —10 e.u. as reported for 
U stacks 0. This choice brackets the range of reported 
parameters and will allow us to set tentative boundaries 
of the experimental conditions to explore. The terminal 
monomers of the domain experience stacking interaction 
with one neighbor rather than two. The reduction in 
their configurational entropy is possibly weaker. To al- 
low for these two effects we assign each terminal monomer 
with an additional free energy Af t . The correspond- 
ing Zimm and Bragg parameters are s = exp(— Af/kT) 
and a = exp(-2A/ t /feT) llj. The 9 vs T melting 
curves are broad leading to 0.5 < a < 1 and suggesting 
weak cooperativity. For simplicity we will assume perfect 
non-cooperative behavior with a — 1. For comparison, 
in helicogenic polypeptides the ith monomer binds the 
monomer i+3 thus giving rise to higher cooperativity sig- 
nalled by much smaller a values, of order 10~ 2 — 10~ 3 [l2fl . 
The distance between the bases in the stacked form varies 
between 0.32 nm and 0.35 nm depending on the mea- 
surement technique and T. In the following we will thus 
assign a value of b — 0.34 nm to the projected length 
of a stacked monomer along the axis of the helical do- 
main. Because of the non-cooperativity of the stacking 
the helical domains are relatively short. While their per- 
sistence length is not known we will assume that it is 
much longer than the typical domain length and thus 
effectively infinite. In contrast we assume that the un- 
stacked domains behave as freely jointed chains and we 
neglect excluded volume effects. Two parameters thus 
characterize the coil-like domains: the effective length of 
an unstacked monomer, a, and the persistence length of 
the coil, A. Neither is well established. A common value 
for a is 0.6 nm |7j while reported values of A, for stacking- 
free chains vary between 0.75 nm and 3.5 nm We 
will utilize A = a and A = 3.5 nm. The free energy of an 
unconstrained homonucleotide within this model is [f|: 

j— = -9\ns + y^2[p s (n)\np s (n) + p u (n) In p u (n)} 

n 

- ^[(^+¥2)P s (n) + W+¥2')P«(n)]- (1) 

n 

The first term allows for the excess free energy of the 
stacked bases. The next two terms specify the mixing 
entropy arising from the polydispersities of the stacked 
and coil-like domains where p s ( u )(n) is the probability of 
stacked (unstacked) domain comprising of n bases. The 
last four terms impose two constraints by use of Lagrange 
multipliers. //* assures the probability normalization 
while H2 imposes the average number of monomers, 
0/y and (1 — 9)/y respectively, in these domains [f|. The 
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FIG. 1: Plots of (R 2 )/Na = 2 l v vs. T for poly(dA) with 
l p given by Eq.Q exhibit a minimum at T m i n = 307 K for 
A = 3.5 nm (full line), T min = 327 K for A = 0.6 nm (dashed 
line) in comparison with (R 2 )/Na vs. T when l p — n/kT (dot- 
dashed line), k was set by equating l p at T m i n for A = 3.5 nm. 
The inset depicts the same plots for poly(rU). 



equilibrium conditions, dFo/d9 = dFo/dy — 0, yield 9 = 
s/(s + l),y = s/(s + l) 2 = 9(1 - 9), p.(n) = (1 - 9)9 n ~ 1 
and p u (n) = 9(1 — 9) n ~ 1 . The average number of bases 
in a stacked and unstacked domain are thus respectively 

(n) s = J2n n Ps( n ) = - 0) and ( n )u = ^nPuW = 
1/9. 

The simplest characteristic of stacking is the T depen- 
dence of the mean square radius of the chain, (R 2 ), as 
determined by the end-to-end distance R. Assuming that 
the rod and coil segments are freely jointed, the con- 
tributions of stacked and unstacked monomers to (R 2 ) 
are independent. The (1 — 9)N unstacked monomers 
in coil domains with a persistence length A constitute 
Na(l — 9)/2X freely jointed segments of length 2A con- 
tributing (2A) 2 iVa(l -0)/2A to (R 2 ). Of the 9N stacked 
monomers, the ones that form domains incorporating n 
bases contribute N9p s (n)/(n) s freely jointed segments of 
length nb. Altogether, 

(R 2 ) = N2\a{l -6) + Nb 2 9 (n 2 ) s /(n) s (2) 

where (n 2 ) s /(n) s = (1 + 9)/(l - 9) 14]. The qualita- 
tively important feature of (R 2 ) is a minimum (Fig.l) at 
#min = 1 — \/2fr 2 / (2Aa + b 2 ). This effect disappears if we 
ignore the differences between the size of the monomer 
in the two states (2Aa = b 2 ). It reflects a competition 
between two contributing processes: (i) the number of 
effective monomers increases with T because the number 
and the size of the helical domains decrease, and (ii) the 
size of the chain at low T is dominated by a single heli- 
cal domains whose length scales with N rather than with 
N 1 / 2 . For the values of a and b we utilize the minimum is 
attained at 9 = 0.769 or s — 3.32 for A = 3.5 nm and at 
9 = 0.474 or s = 0.9 for A = 0.6 nm. Thus for poly(dA) 
the corresponding T m - ln values are T m ; n = 33.6°C and 
T min = 53.5°C while for poly(rU) T min = -55°C and 
T m i n = 2.5°C. A closely related effect occurs in helico- 
genic polypeptides as they undergo a cooperative helix- 
coil transition |l2j. 
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ssNA chains are occasionally considered as semiflexi- 
ble, worm-like chains. In this case, it is implicitly as- 
sumed that the properties of the chain along its back- 
bone are uniform. In the limit of iV > 1 scmiflcxi- 
ble chains obey (R 2 ) — 2l p L where l p is the persistence 
length and L = Na is the contour length. Within this 
model l p — n/kT where k = const is the bending modu- 
lus of the chain. If the behavior of the ssDNA is analyzed 
within this framework while the chain obeys the stacking 
model, the T dependence of l p thus extracted is given by 

2l p = 2A(1 - 6) + (b 2 /a) 9(1 + 9)/{l - 9). (3) 

At high T, when 9 — > 0, l p « A while for low T, when 9 —> 
1 it diverges as l p ~ b 2 /a(l - 9) ~ b 2 cxp(-A//fcT)/a. 
This T dependence is markedly different from the 1/T 
behavior predicted by the worm-like chain model (Fig.l). 

A related signature involves cyclization reactions. Cy- 
clization reactions require the two ends of the chain 
to be within a certain capture radius, r c <C (R 2 ) 1 ' 2 . 
The thermodynamics of the ring formation are de- 
termined by P(R)dR, the probability for the end-to- 
end distance of the chain to be in the range R to 
R + dR. When N > 1, P(R) of flexible homopoly- 
mers, behaving as freely jointed chains with constant 
monomer size, assumes a Gaussian form 15]: P(R) — 
AttR 2 (2tt(R 2 (9)) /3)- 3 / 2 exp(-3i? 2 /2(R 2 {9))). This re- 
sult applies to ssDNA in the limit of 9 — » 0, when the 
effect of the stacking is negligible. Clearly, it is wrong 
when 9 — > 1 and the configurations are dominated by 
a single, long stacked domain. For 9 > 0, the Gaus- 
sian form is valid provided N 3> 1 and yN 3> 1 i.e., 
the stacked monomers form a large number of domains. 
The polymer may then be considered as a freely-jointed 
chain whose effective monomers are rod-coil diblocks of 
varying sizes. The Gaussian form is applicable in this 
regime since the probability distribution of lengths of the 
rod-coil "monomers" , while unspecified, is identical for 
all rod-coil diblocks [16J. The cylization equilibrium is 
ruled by the elastic free energy arising from constrain- 
ing R, F el = -kTlnP(R) = -kT\n{R 2 /{R 2 (9))V 2 ) + 
2>kTR 2 /2(R 2 (9)) + cste. This F ei has negligible effect 
on 9 and y because its contribution to the equilibrium 
conditions arises from R 2 /{R 2 (9)). For R < N 1 / 2 a the 
corresponding terms scale as 1/N and are thus negligi- 
ble. To obtain the precise cyclization penalty it is im- 
portant to allow for the weighted contributions of all the 
configurations with R < r c . Since for r c <C {R 2 {9)) 1 / 2 
the exponent in P(R) is of order unity, the fraction of 
cylizable states within the freely jointed chain model 
is J^P(R)dR ~ [r c /(R 2 (9)) 1 / 2 } 3 . For self avoiding 
chains (R 2 (9)) 1 / 2 is replaced by the 9 dependent Flory 
radius 17]. Altogether the cyclization entropy for r c <C 
N x l 2 a is S cyc (9) = ik\n[r c /(R 2 {9)) 1 ' 2 ] and the equi- 
librium constant for the cyclization reaction is speci- 
fied by kT\nK cyc = AH - TS cyc {9) where AH is the 
binding enthalpy of the terminal groups. The activation 
free energy for cyclization, AF^ may be identified with 
-TS cyc {9) ~ kT\n(R 2 (9)) 3 / 2 3]. This AF* exhibits a 



minimum at m i n and thus at the corresponding T m ; n . We 
should emphasize that this analysis is not valid when the 
chain contains a large stacked domain comprising most of 
the monomers. In this case the bending of the chain can 
induce melting of the stacked domain thus introducing a 
coupling of 9 and the cyclization: The elastic free energy 
of fully stacked chain of length L = l p = n/kT forming 
a ring of radius R = L/2tt is kL/2R 2 = 2n 2 kT w 20/cT 
while the reported stacking free energy per base at 25°C 
is in the range of 2 to 8kT. This rough argument suggests 
that the chain can lower its free energy by "melting" a 
few stacks thus avoiding the bending penalty. 

When the reaction between the terminal groups is dif- 
fusion controlled the cyclization rate constant assumes 
the form k cyc ~ 1/t where r is the longest character- 
istic time of the chain |18| . If hydrodynamic interac- 
tions are neglected r = where tr w (r] s a 3 /kT)N 2 is 
the Rouse time and n s is the solvent viscosity. Allow- 
ing for hydrodynamic interaction leads to r = tz where 
t~z ~ rjsR 3 /kT is the Zimm time. When excluded volume 
interaction are negligible tz ~ A 3 / 2 while in the opposite 
case Tz ~ N 3v where v = 0.588. Experimental studies 
of cyclization of synthetic polymers in non-aqueous solu- 
tions are consistent with k cyc ~ N~ 3 / 2 . Our discussion 
suggests thus that the cyclization rate constant of long 
ssDNA will follow k cyc - (R 2 {9))- 3 / 2 with a maximum 

at $min- 

The T dependence of R, K cyc and k cyc is recovered if 
we consider the ssDNA as a worm like chain with a pre- 
sistence length l p given by 0. This picture fails qual- 
itatively when considering the extension force law. To 
obtain it, we augment Fq with the elastic free energy, 
F e i, of a freely-jointed chain subject to tension / 0|. A 
Kuhn length 2A is assigned to the coil domains while each 
of the rigid helical domains orients as an independent ef- 
fective monomer. Allowing for the size distribution of 
the stacked, helical domains we obtain [j| 

= - ^—r^-£int (Sx) - y Ps (n)Ci n t {inx) (4) 

n 

where 5 — 2X/a, 7 = b/a, Ci„t(x) — ln[sinh(x)/a;] and 
x = fa/kT. Minimization of F = Fo+F e i in the constant 
/ ensemble yields 6, y, p s ( u )(n) [2£j and the force law in 
the form 

^- = (1 - 9)C{5x) + 7 coth( 7 x) - - ^. (5) 

Na sinh(7a;j x 

At low T the force laws obtained from this model exhibit 
a smoothed plateau associated with the enhancement of 
stacking by the applied force (Fig. 2). Initially, the ex- 
tension lowers the entropy of the chain thus favoring the 
ordered, stacked domains. Eventually, further stretch- 
ing enforces melting of the stacks in order to release the 
stored length (b = 0.34 nm vs a = 0.6 nm). For our choice 
of parameters, poly(dA) exhibits a pronounced plateau at 
T = 25°C while for poly(rU) significant deviations from 
the freely jointed chain appear below T — 10°C. No de- 
viations are expected for poly(dT). Modelling the chain 
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FIG. 2: Plots of / vs R/Na for poly(dA) at T = 25° C as 
obtained from the stacking model Eq.@ with A = 3.5 nm ( — 
-), from the freely jointed chain model with l p given by Eq.© 
(xxxx) and with l v = 3.5 nm (oooo). Insets: (a) depicts the 
same plots for poly(rU) at 10°C. (b) Schematic picture of 
ssDNA with stacked blocks. 



as a freely jointed chain with a T dependent Kuhn length 
© does not recover this plateau because it does not allow 
for the coupling between the stacking and the extension. 
Accordingly, single molecule stretching experiments pro- 
vide a stringent test allowing to distinguish between the 



conflicting views on the effect of stacking on the elasticity 
of ssDNA. 

Straightforward modification of models for the helix- 
coil transition in polypeptides allowed us to study the ef- 
fects of stacking interactions on the configurational and 
elastic properties of homonucleotides. In particular, we 
investigated the dependence of (i? 2 ) 1 / 2 , K cyc , k cyc and 
the force laws, / vs R, on T. The most noticeable effects 
are extrema in the plots of (R 2 ) 1 ' 2 and k cyc vs. T as well 
as the appearance of a plateau in the force law. The nu- 
merical results are based on reported values of Ah, As, 
a, b and A as obtained from experiments. In confronting 
our predictions with future experiments it is important 
to note that the reported parameters vary widely with 
the measurement technique and the experimental condi- 
tions. They also depend on the model used to analyze the 
data. With this caveat in mind, the qualitative effects we 
discuss provide powerful diagnostics for the coupling of 
stacking interactions with the elastic properties of ssNA. 
These predictions are meaningful because spectroscopic 
evidence indicates significant stacking in poly (d A) and 
poly(rA) at 20°C, irrespective of the precise values of 
Ah, As, a and A. The results presented above can also 
be used to test the performance of various sets of Ah, 
As, a, b and A in recovering the observable (R 2 ) 1 / 2 vs T 
or / vs R plots. 
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